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I. INTRODUCTION 

Now properties of multiparticle entangled states, their preparation and application are the subject of extensive 
discussion. The desired state of a physical system can be prepared either by projective measurement or as a result 
of evolution. For atomic systems, both methods have already been implemented in experiments. In particular, 
two atomic ensembles were used in |lj to create an EPR pair by projective measurement. The latter method was 
demonstrated in several studies: entangled states of alkali ions were generated via Coulomb interaction Q, neutral 
Rydberg atom were used to create an EPR pair in a microLaser setup 1, and resonant dipole - dipole Interaction 
was used for entangling neutral atoms in an optical lattice The most popular methods for preparing entangled 
photon states are still mostly based on parametric down-conversion. For example, an entangled state equivalent to a 
three-state quantum system (qutrit) was prepared and examined by using quantum state tomography in [5|. These 
examples suggest that an entangled state of two systems can be prepared experimentally by using a certain interaction. 
Systems of this kind are well studied. With regard to applications, it is important to know how entanglement can be 
utilized and its robustness against decoherence. In this respect, of special interest are multiparticle systems, whose 
entangled states are characterized by much more complicated and diverse behavior. 

Previous efforts were mainly focused on analysis of entanglement between several particles. In particular, the W class 
of tripartite entanglement defined in [fj includes the symmetric three-photon polarization entangled state implemented 
in the experiment reported in Q ■ An extension to four qubits was proposed in Q , where nine inequivalent classes were 
distinguished that cannot be connected by local operations and quantum communication. Studies of multiparticle 
systems are relatively few, being focused on entanglement criteria and application to problems in quantum information 
theory. Whereas the Peres - Horodecki criterion for bipartite entanglement found in \$ was applied to a real physical 
system in no operational criterion is known for ent angl ement in the general case, and various approaches are 
often used. In the concept of entanglement molecules 0] was introduced to propose a classification using graphs, 
in which particles and classical or quantum correlations represented, respectively, by vertices and edges connecting 
pairs of vertices. Graphs of this kind can be used to describe both pure and mixed entangled states and distinguish 
several classes differing by topological properties of the graphs. In [Rj, symmetric states (including Dicke states) were 
studied by using several entanglement measures (entanglement entropy, negativity, and entanglement of formation) 
defined by the eigenvalues of a partial transpose of the density matrix. A numerical analysis was performed to find 
that symmetric states are robust to particle loss even if the number of particles is large (up to 10 3 ). Note that 
the calculation of eigenvalues is a difficult task, because the dimension of an ensembles Hilbert space exponentially 
increases with the number of constituent particles. Owing to their robustness, symmetric states can be used in such 
applications as cloning and telecloning protocols for quantum information transmission , quantum key distribution 
, and quantum teleportation or dense coding jT?3 | . The formulation of two models of a one-way quantum computer 
using measurements on multiparticle entangled states 0, has strongly stimulated studies of the properties of 
multiparticle systems, in particular, Ising- and Bose - Hubbard-like models. 

The present study focuses on the Dicke states arising as a result of collective interaction of many atoms with 
electromagnetic field ^J, which has been analyzed in numerous studies (e.g., see 113] )■ This system exhibits many 
physical properties of interest for quantum information processing. Photon trapping in chain configurations of atoms 
was considered in |2l( . When the system is placed in a cavity, this effect reduces the photon escape rate and increases 



the decoherence time of the cavity mode. In J22(, this effect was used for generating W states and anticloning 23], 
which can be implemented with high fidelity by means of photon trapping. In those studies, only single-photon traps 
and single-photon initial states were analyzed. Here, we consider the more general case of multiphoton processes, 
assuming that the photon statistics is arbitrary. 

The main questions addressed below arc the following: what types of entangled states are produced by interaction 
between atoms and field? What states can be prepared from independent atomic ensembles entangled with a photon? 
How can these states be utilized? We consider resonant interaction between narrow-band light and an ensemble of 
identical two-level atoms coupled to a common heat bath. The analysis is restricted to a simple model of radiative 
decay. Multiphoton processes, such as Raman scattering, are described in terms of effective Hamiltonians, which can 
be obtained by unitary transformation |21| . The behavior of an atomic system interacting with light characterized by 
arbitrary photon statistics is analyzed by using perturbation theory in the interaction strength for arbitrary statistics 
oh light particularly for Gaussian, coherent, and squeezed states. We find that weak single- mode squeezed light is 
required to create multiparticle entanglement between atoms. As distinct to the case considered in |25j, the steady 
state discussed here is robust against atomic decay. When decay is neglected and analysis is restricted to a single- 
photon initial state, simple exact solutions describing exchange of excitation between the field mode and atoms can be 
obtained (2(| . These solutions can be used for generating and transforming symmetric Dicke states and for processing 
and storing quantum information. The optical schemes for projective measurement considered here can be used to 
generate entangled states of atomic ensembles. An EPR entangled pair of macroscopic ensembles was created in an 
experiment 0. The new states produced in our schemes have hierarchical structure, thus differing from the cluster 
states introduced in [^J as a resource for one-way computing. 

The paper is organized as follows. First, we formulate a basic model and write out the second-order perturbation 
solutions obtained by taking into account radiative decay. These solutions are then used to analyze the states of the 
atomic system corresponding to various photon statistics. Exact solutions obtained under certain initial conditions 
by neglecting radiative decay are used to describe generation and transformation of symmetric Dicke states. Finally, 
we consider optical schemes for preparing entangled states of atomic ensembles by projective measurement. 



II. BASIC EQUATIONS 

In the dipole approximation, the ensemble of N identical, but distinguishable, two-level atoms interacting with 
electromagnetic field is described by the Hamiltonian 

h = m, 

■d = ^gka-kSl - h.c, 

k 

where gu = (fiw / t/2eoi 3 ) 1 ^ 2 (/i, e^) is the coupling constant, /i is the dipole transition matrix element, is the 
polarization vector for the mode with wave vector k, and a\, are photon creation and annihilation operators, 
St = ^2 a sio(a) exp(ifcr a ) is the atomic operator for the atom located at a point r a (x,y = 0, 1, where and 1 denote 
the ground and excited states, respectively). When analysis is restricted to interaction with a single resonant mode, 
Sk can be replaced with Sk=o, which makes it possible to treat an atomic ensemble occupying a spatial region as a 
point like object. Then 

# = S W B-S 01 B\ (1) 

where Siq = |l)a(0|, B = ga. Effective Hamiltonian JQ) is used here to describe not only interaction with a single 
resonant mode, but also multiphoton processes, such as Raman scattering. In the latter case, we set B = /a^a^, and 
assume that the photon frequencies loa and los satisfy the relation u> — lja — u>s_j where lo is the atomic transition 
frequency. Hamiltonians of this form can be obtained by unitary transformation [24J . 
The density matrix p of N atoms interacting with field obeys the master equation 

^- t p=[#,p]+£p, (2) 
where relaxation is represented by the Lindblad superoperator 



a 



C a = — j[soi(a)s w (a)p - s 10 (a)ps 01 (a)} - y[sio(a)soi(a)/9 - s i{a)ps 10 {a)} + h.c. (3) 

This representation corresponds to the model of purely radiative decay with longitudinal and transverse decay rates 
7 = 7j_ + 7f and j± , which satisfy the relation 7^ = 7/2. In general case j± > 7/2 since 7^ should be replaced by 
7^ + K, where k is a dephasing collision rate. 

Effective Hamiltonian JQ) may involve many field modes with u>k differing from the atomic transition frequency by 
5uk and occupying a frequency band of width Aw. If Aw, 5u>k <C 7_l, then we can consider a narrow- band radiation 
field and make use of resonance approximation. Otherwise, the field must be described in terms of multiple-time 
correlation functions. Solution of Eq. |2"|l is a difficult task. To describe the interaction between single atom and field, 
the following equation for the density matrix p a = Tr' a p is derived from J5J by tracing over all atoms except for one: 

d 

-7T,Pa = [#a,Pa] + ^aPa + N(N - l)Sp a > O ', Paa' , (4) 

at 

where i? a = s\o[a)B — h.c. and p aa > = Sp' aa ,p is a two-particle density matrix. The right-hand side of contains 
a multiparticle contribution proportional to N(N ~ I), because the density matrix p aa / does not commute with the 
field operators. This leads to the Bogolyubov- Born-Green-Kirkwood-Yvon chain of equations for the multiparticle 
density matrices p a , p aa ' , Paa' a" , ■ ■ ■ • In physical terms, this means that fluctuations of quantized electromagnetic field 
induce correlation between atoms. If the field is assumed to be classical and noise- free, for example, a coherent state is 
considered, then the interaction will not give rise to any correlation, and the initially uncorrelated atoms will remain 
mutually independent. In what follows, we use (J2J to analyze interactions that can be used to generate symmetric 
Dicke states. 



III. DICKE STATES 



First, we define symmetric Dicke states and introduce a representation of symmetric Dicke states that demonstrates 
their relation to the collective interaction processes. The Dicke states are eigenstates of the operators J z and J 2 = 

J z \jma) = m\jma) 
J 2 \jma) = j(j + l)\jma), (5) 

where [J s , J p ] = ie sp dJd- For example, operators J Sl s — x,y, z can be represented by Pauli matrices J s — (1/2) J^ fc a s k, 
&sk,s — x,y,z. Indexes j and m are integer or half-integer numbers |m| < j, maxj = N/2. If j = N/2, then the 
states are symmetric, and the quantum number a introduced to lift degeneracy can be omitted. For h excited atoms 
h = m + N/2, the states can be represented as 

\j = N/2,m) = \h;N)=Y l P^i,h,...,lh,0h+i,...,0), (6) 

z 

where P z is one of the = N\/(h\(N — h)l) distinguishable permutations of particles. 

The vector \h; N) describes an atomic of h excited atoms and it is normalized as (h; N\h; N) = . Symmetric states 
of a multiparticle system arise when interaction is described by collective operators of the form S\q — 1 1 ) « (0 1 . 

\h; N) = (1/W)S£,|0;JV>. (7) 

If ft, = 1, then one finds that 

|1; N) = |10. . .0) H + |00...1). (8) 

Since the wavefunction \h;N) is not factorizable, it represents an entangled state. In terms of correlation between 
particles, it is substantially different from other entangled states. For example, in the Greenberger HorneZeilinger 
(GHZ) state GHZ = (1/2)(\0)® N + \1)® N ), the correlation of any M particles (M < N) is classical. In particular, the 
density matrix corresponding to the state |1; iV)(l; N\ of a group of M particles is p(M < N) = iV -1 |l; M)(l; M\ + 
(N— M )N~ 1 10; N) (0; N\. The corresponding von Neumann entropy depends on the relative particle number p = M/N: 
S(p(M < N)) — —p\ogp — (1 — p) log(l — p). When p = 1/2 the entropy achieves its maximum 1. If M = 2 we 
can apply the necessary and sufficient separability criterion proposed in According to this criterion, the state is 
inseparable (entangled) if the density matrix partially transposed over the one of the atoms has at least one negative 



eigenvalue. In the case considered here, one of the four eigenvalues {1/A; 1/A; (A — 2)(2A)~ 1 [1 ± ^/l+AJ (A — 2) 2 ]} 
is negative. Note that the behavior of correlation between M particles depends on p = M/N . As the total particle 
number A increases, p — > and the correlation vanishes, since their state becomes pure as p(M < A) — > |0; N)(Q; N\ 
In what follows, we make use of the following equalities: 

S 01 |0;A) =0, 
S 10 \h;N) = (h + l)\h+l;N), 
S 01 \h;N) = (iV-/i + l)|/i-l;AT), 
5 01 S 10 |ft;iV) = (ft+l)(JV-/i)|fc;JV), 
SioS , oi|/i;AT) = ft(iV-/i + l)|/i;iV). (9) 



IV. SECOND-ORDER PERTURBATION THEORY 



To solve Eq. J5J), we use perturbation theory in the interaction strength: 

p = p^+p^+ p W + ..., (10) 

Here, the zeroth-order approximation p(°) is the steady-state solution of with i? = 0: p^ - 1 = |0)(0| ® p/, where the 
density matrix pf represents the modes and |0) = |0)® corresponds to the ground state of all atoms. The operators 
p( k \ k = 1, ... satisfy the equations 

= (11) 

subject to the initial conditions p^ (0) = 0. 

The analysis that follows is restricted to second-order perturbation theory, which is sufficient to obtain statistical 
characteristics of the excitation field. The matrix equation for p( 2 ) is 

(l fc ,l m ;A|-|p( 2 )|0;A) = -2 1± (l k ,l m ;N\pW\0) + (l k ,l m ;N\R\0) 
at 

(l fc ; A||-p( 2 )|l m ; A) = -2 7± (l fc ; A|p( 2 >|l m ; A) + (l fc ; N\R\l m ; A), k + m 
at 

(l fc ; 7V|^p (2) |l fe ; AO = - 7 (l fc ;iV|p( 2 )|l fe ; AO + (U; N\R\l k ; A) 
(0; N\§-/ 2) |0; A) = 7 ^(l fe ; A|p( 2 ) |l fe ; TV) + (0; A|i?|0; A). (12) 

where sio(fc)|0; A) = |lfe;A) and sio(fe)sio(p)|0; A) = |l/.,l p ;A) represent the states in which only the kth atom 
is excited and only the kth and pth atoms are excited, respectively. The nonzero matrix elements of the operator 
R = [#, pW] are 

(l k ,l m ;N\R\0;N}=2K{t)B 2 p f , 
(0; N\R\0; A) = -n(t)N{B^ Bp f + pfB^B), 

(l fc ; A|i?|l m ; A) = 2 K {t)Bp sf B\ (13) 

where n(t) = (1/7j_) (1 — exp(— 7 j_t)). For purely radiative decay, 7 j_ = 7/2 and the second-order perturbation theory 
yields 

p = |0) (0| ® p f + k[\1\ A) (0; A| ® Bp s/ + ft.c] + k 2 [|2; A) (0; A| (g) B 2 p 7 + h.c] 
-A 7 /C|0; A)(0; A| [fl t Sp / - BpfB^ + h.c] - (1/2)Ak 2 |0; A)(0; A| [B^/O/ + /i.e.] 

+K 2 \l:N)(l;N\(g> Bp f B\ (14) 

where 

,-1 J «.-2r • 71 



£ = 7I 1 <j 7~ 2 [7* + 1 - exp(-Tt)] - — 



This expression is valid to second order if the field is relatively weak: 



Nn 2 (B^B) < 1. (15) 

In the case of interaction with a single resonant cavity mode, we have B = ga and n 2 (B^B) — n/n s , where n s — 
(7J-/.9) 2 is a saturation parameter and n = (a^a) is the mean photon number. Then, (|15f) reduces to the standard 
condition imposed in the case of resonant coupling between the field and two- level atoms: Nn/n s -C 1. Solution 1|14|) 
describes the joint evolution of the atomic ensemble and field starting from an ensemble of ground-state atoms and 
an arbitrary state of the field. 

V. MIXED NONSEPARABLE ATOMIC STATES 

Second-order perturbation theory predicts correlation between atoms depending on photon statistics, i.e., provides 
a framework for describing entangled (inseparable) atomic states. To analyze the properties of the atomic system, we 
use second-order perturbation theory to find the density matrix for a group of M < N atoms, pa(M < N), obtained 
by taking the trace of l|14(l over the field states represented by pf and over N — M particles. The result has the form 

p A {M <N) = |0)(0|[1 - M n 2 (B^ B)] + k[(B) |1; M) (0| + h.c] + k 2 [(B 2 ) |2; M) (0| + h.c] 

W{B^B)\1;M){1;M\. (16) 

Note that the density matrix pa(M < N) describes a mixed state of the atomic ensemble. Unlike the density matrices 
for symmetric Dicke states Pa(M < N) is independent of both N and p = M/N. Therefore, the correlations 
between M < N atoms are identical and are independent of the total particle number. This implies that the state is 
robust to particle loss. 

The atomic density matrix cannot be factorized because of the correlation depending on photon statistics. Consider 
two atoms described in terms of their respective observables c\ and c 2 such that [ci; c 2 ] = 0. Setting M = 2 in 116|l . 
we have the two-atom density matrix 

p A {2) = |00)(00|(1 - 2k 2 (B^B)) + /t(B)(jlO)(00| + |01) <00j + /i.e.) + k 2 (B 2 ) (|11)(00| + h.c)j 

+ K 2 (S t S)(|10> + |01>)((10| + (01|). (17) 

Using l|17f) we find that the covariance of the operators ci,C2 is determined by the electromagnetic field variance: 

(cic 2 ) - ( Cl )(c 2 ) - n 2 [((B 2 ) (B) 2 )(0| Cl |l)(0|c 2 |l) + ((B^B) (B^(B))(l\ Cl \Q)(Q\c 2 \l) + c.c]. (18) 

If the field is not fluctuating in the sense that its variances are zero, i.e., (B 2 ) — (B) 2 = etc. (which is true in the 
present case, e.g., for a coherent state), then there is no correlation between atoms. Suppose that Ck(k = 1,2) are 
dipole operators: Ck — — p(soi(k) + sia(k)), where the matrix element p is real. Then the correlation between 
two dipole moments depends on photon statistics. We define the quadrature operator X f = B^ exp(i6) + h.c. Then 
()18(l implies that the covariance of the dipole moments is determined by the variance of the quadrature operator 
normally ordered with respect to the field operators B and B^ at 9 = 0: ((iiffe) — (^i)(^2) = P 2 k 2 Dn, where 
£>jv = (Xf) — {Xf) 2 — ([B, B^]). For coherent states, the variance is Dn = 0. The dipole moments are correlated 
both for a squeezed-state field (with D n < 0) and for field in a classical state (with Dm > 0). 

The necessary and sufficient condition for inseparability of a mixed state is provided by the Peres Horodecki 
criterion 0, which is valid for systems with Hilbert spaces of dimension 2x2 and 2x3. In the case considered here, 
the state of a two-atom system described by pa (2) is inseparable (entangled) if at least one eigenvalue of the density 
matrix partially transposed over the variable of atom 1 pa(2) Ti is negative. As example, we consider light in Gaussian 
and squeezed states. 

For a Gaussian field ((B) = (B 2 ) = 0, expression l|17fl reduces to the density matrix describing a superposition of 
the ground and mixed states: p A (2) = a|00)(00| +&[(|01) + ]10>)«01| + (10|)], where a + 2b = 1 and a = 1-2k 2 (B^B). 
The eigenvalues of pa(2) Ti are 

Since ^Ja 2 /A + b 2 a/2, in the approximation considered here, we have the eigenvalues: {b, b, a,0} i.e., a separable 
state. 



Consider the case of resonant interaction with single- mode squeezed light (B = ga) generated, for example, by a 
parametric oscillator. A simple model of the oscillator is defined by the effective Hamiltonian H — ih(f /2)(a^ 2 — h.c). 
The solution is a — aocoshr + sinhr, where r = fr is the squeezing parameter, r is the normalized length 
of the nonlinear medium, and ao,aJ and denote the input field operators. For the initial vacuum state, (a) = 0, 
(a 2 ) = (a^ 2 ) = coshr sinhr, (a^a) = sinh 2 r. In this case ()17f) reduces to the following the two-atom density matrix 



p A (2) = 1 00) (00 1 [1 — 2k 2 (B^ B)} + k 2 ^(B 2 } (J 11) (00 1 + |00)(ll|j + h.c. 

+ K 2 (B t B)(|10)(10| + |01)(10| + |10)(01| + |01)(01|). (19) 

The four eigenvalues of p 7 ^ (2) are 

A = I 0; 1 - — sinh 2 r; ± — cxp(±r) sinh r \ . (20) 
L n s n s J 

To be specific, we set r > 0, i.e., consider the state squeezed with respect to canonical momentum or phase. In this 
case, (— l/n s ) sinhr exp(—2r) < 0. However, it is clear that the degree of squeezing is low, because the approximations 
used here imply that 



sinh 2 r 



« 1. (21) 



Thus, the state of the atomic system is inseparable. This behavior is explained as follows. Fluctuations of light give 
rise to correlation between atoms, which leads to two-atom coherence. When condition 1|21|) holds, this coherence 
plays the key role. Since absorption is weak, the system is almost entirely in the ground state. As distinct to the case 
of Gaussian statistics, the density matrix has the form pa(2) ~ 1 00) (00 1 + k 2 [(B 2 )|11)(00| + h.c.]. 

Note that the following two observations can be inferred from this example. First, a steady entangled atomic state 
can be created by using weak squeezed light, which looks promising from an experimental perspective. Second, the 
entire ensemble cannot be interpreted as separable, because any pair in a group of M < N atoms is entangled, i.e., 
the quantum correlation of the ensemble as a whole is robust to particle loss 

Since no reliable universally applicable criterion is known for multiparticle entanglement, we apply the Peres 
Horodecki criterion to two two-level subsystems and found that any pair of atoms in the ensemble can be insep- 
arable, which gives reason to interpret the state of the entire system as inseparable. 

Note also that spurious entanglement may be predicted by perturbation theory |28j|. In that study, an example of 
expansion of the product of two wave functions in terms of a common classical parameter was considered in which 
individual summands represent entangled states. However, if entanglement entropy is used as a measure, then we have 
initially independent systems, because the entropy is cither quadratic in the small parameter or zero in arbitrary-order 
perturbation theory. Note that physical implementation of such entangled states, i.e., preparation of an independent 
state of a pair of entangled particles, requires projective measurement in an entangled basis. The present analysis 
also relies on perturbation theory, but we deal with a different situation in both physical and formal sense, in 
which interaction between particles gives rise to correlation. The wavefunction obtained in first-order perturbation 
theory is not factorizable, and the corresponding entanglement entropy is zero to the corresponding accuracy. This 
result is physically plausible, because there is no correlation in the first-order perturbation theory. In our analysis, 
entanglement is predicted by second-order perturbation theory, which describes real emission and absorption processes 
result in correlation. In this order of perturbation theory, the existence of quantum correlation is substantiated by 
entanglement criteria consistent with approximation accuracy. 

VI. EXACT SOLUTIONS 

Radiative decay can be neglected in when evolution over a time t <C 7" 1 is considered, and the behavior of the 
entire system is described by the wavefunction <j>{t) = exp(—ih~ 1 Ht)((j)A <8> 0/), where the initial states of the atoms 
and field are assumed to be uncorrelated. Then, simple solutions can be obtained under certain initial conditions. 

Consider the mixing of modes a and b described by 



H = ihf{a i bS — ab^S^), 



(22) 



where S = SicS^ = Sqi. If analysis is restricted to single-photon Fock states of the modes <f>f = c|01) a & + e|10) Q & 
exact solutions can be written as 

eBcp{-»» _1 frt}(c|01) O 6 + e|10) o6 ) ® $ A = c{|01) cos[tfVsS^} + \io)S^-^= am[tfVsS^}\ O cj) A 

+e{ - |01>S'-^=sin[i/V^t5] + |10) cos[i/\/si\S] } ® <f> A . (23) 

In the case of a single-photon process described by the Hamiltonian 

H = ihg{aS~a ] S ] ) (24) 

there also exist simple solutions. For example, 

exp{-ift -1 frt}(c|l) <g> |0; JV) + e|0) <g) |1; JV)) 
= c{cos[ 5 /ViV]|l) <g> |0;JV) + sin[ 9 /ViV]|0) ® |l;iV)} (25) 

+e{ - V7Vsin[ ff /ViV]|l) <g> |0;JV) + cos[<?/ViV] |0) <8> |1;JV)}, 

where \h; N) = 10)®^, h = 0, 1 represents the ground state of the atomic ensemble and a symmetric Dicke state defined 
in accordance with (jHJl. These solutions are valid only under the restrictions imposed above on the initial states. They 
describe exchange of excitation between the cavity mode and the atoms. 

VII. GENERATION AND TRANSFORMATION OF SYMMETRIC STATES 

Now, we use the exact solutions written out above to analyze the evolution of symmetric Dicke states \h;N) in 
single- photon and wave- mixing processes. 

First, consider the case when the spatial inhomogeneity of the field within the region occupied by the atomic 
ensemble can be neglected. Setting (|2*5|) 4>a — \h;N), we use © to obtain 

(a\01) +/3|10)) <8> \h;N) a{ cos 6^01) <g> \h;N) + JlEEL S m6 h \10) <8> \h+l;N)} 

+l3 { ~ \/ N ~h +1 sin6) ^ 01 ^ \h-l;N) + co8 6' h \10) ® \h;N)} (26) 

where Oh — tf-J (h + l)(N — h), 6' h — tfy/h(N — h + 1). Relation i|26|) entails possibilities of preparation of an 
entangled from ground-state atoms |0;iV) — > |1;7V), and transformation of entangled states by changing the number 
of excited atoms \h; N) — ► \h ± 1; N), including disentanglement: \h; N) — > \h — 1; N) — ► . . . |0; AT). 

Note that exact solutions (|25|l and (|26|l describe state swapping, which can be used to map the state of light onto 
atoms in order to store it in a long-lived atomic ensemble, i.e., to implement quantum memory. In particular, an 
unknown superposition of photons can be transferred to atoms and back by using the following transformation entailed 
by 

(a\l) + m) ® |0; N) ±3 |0) ® (<*-L\l; N) + (3\0; N)) . (27) 



Solutions (|25J) and (|26|l make it possible to take into account the spatial configuration of atoms in the ensemble. For 
example, consider the interaction between a one-dimensional array of atoms located at points x\, . . . , Xn and a single 
photon described by Hamiltonian l|24|) with S — s io(p) exp[ifcaj p ], where s±o(p) — |l) p (0| corresponds to the atom 
located at x p , p = 1, . . . , N. Using (J2EJl, we can show that 

|1)®|0;JV) ->cos6>|l)® |0; N) + sin0|O> <g> r?jv(l), (28) 

6 = tgVN, 

T) N = (l/VN)\e tkxi \10...0) + ...e ikXN \0...01)}. (29) 

Expression (|29|l implies that an array of entangled atoms is created when 9 = tt/2. Note that r)N is the Dicke state 
with j = m = N/2 — 1 only if ^ exp[ikx p ] = 0. 



VIII. ENTANGLED ATOMIC ENSEMBLES 



Solutions l|23l) . I|25|l) imply that a photon and an atomic ensemble are entangled via interaction. If photons are 
entangled (e.g., by projective measurement) in a combination of such independent systems, then the atomic ensembles 
will become entangled. We consider optical measurement schemes based on this method, known as entanglement 
swapping. The key resources used in these schemes are set of atomic ensembles correlated with respective photons, 
beamsplitters, and single-photon detectors. The analysis that follows is restricted to schemes in which only specific 
single-photon output is recorded. 

As an initial state, we use the EPR pair 



Z(W) = a\0) 



f 



■b\l)f®\W), 



(30) 



where Fock states are denoted by the subscript f, \W) — |1; N)/y/N, |0) = |0; N). It is generated by the mode mixing 
described by (|22[1 . where the mode b is a classical wave. The state of n independent identical ensembles entangled 
with respective photons is represented by the product 



Z n {W) = Z(Wf 



1 10 . . . 0)/ ® \W0 . . . 0) + . . . 1 00 . . . 1)/ <g> 1 00 . . . W) 



(31) 



As illustrated by the figure, the photons associated with atomic ensembles are injected into a system of n — 1 
beamsplitters with n input ports and n output ports. Each beamsplitter performs the transformation |01)/ — > 





FIG. 1: (a) Scheme for generating entangled states of atomic ensembles, (b) Preparation of entangled states by correlation of 
photocounts recorded by two schemes. 



Cfe|01) / + Sfc 1 10) /, 1 10)/ — > — Sfc|01) / + Cfc|10) /, where c\ + s\ = 1, k = 1 ... n — 1. The scheme is described by a 
unitary operator U n f and characterized by the following property. There exist an input port optically coupled to 
every output port and an output port optically coupled to every input port. In Fig. la, the latter is output port 1. 
We call it the optical output port, and the corresponding detector is called the output detector. The scheme performs 
the transformation 

[/„/|1...0)/ = t 1 |1...0)/ + --- + t n |0...1) / , 

U~}\1 . . .0)/ = n |l . . .0)/ + • • • +r n |0. . . 1)/, (32) 

where the coefficients tk, Tk, k = 1, . . . n are determined by the transmittances and reflectances of the beamsplitters, 
andE fe ^ = E fc T fe 2 = l. 

If output detector detects a photon which corresponds to the state |1/) = |1 . . .0), then with the probability 

Prob(l) = la^&l 2 (33) 
entangled state of atomic ensembles will be prepared 



(l f \U nf Z n (W)/y/Prob(l) = T) n (W), 

r)n(W) = qi \W...0)+q n \0...W), 



(34) 



This scheme has the following property. Since the coefficients qi, ..q n are completely determined by the transmittances 
and reflectances of the beamsplitters, weakly entangled states Z(W) can be used to prepare highly entangled states 
atomic ensembles. 

Let us consider several particular cases. If n = 2, then qi = c%,q2 = si, and we have an EPR pair of the 
form 772 (W) = EPR{W) = a\W0) + s^OW). When c x = s x = l/\f2 it is maximally entangled. If n = 3 and 
qi = C1C2, q2 = — S1C2, 93 = S2, then we have a W state. If c\ = —s\ = l/V%, C2 — -\/2/3 and C2 = -\/2/3, then 

= = (1/a/3)(|V700) + |0WD) + |0fW)). (35) 

In particular, one can prepare the asymmetric state W(W) = (l/v^|WO0) + (l/2)|0W0) + (l/2)|00W). When N = 1, 
it is unitary equivalent to the GHZ state and can be used as a quantum channel for teleportation or dense coding 

Using correlation between photocounts in a combination of schemes considered above, mixed states of atomic 
ensembles can be prepared, including inseparable ones. For example, consider two independent identical schemes 
S2(X) combined as shown in Fig. lb, with three single-photon detectors in each scheme. If a photon is detected by 
either scheme, then we have the pair of states (lf\S 2 (X) ® (0f\S 2 (X)w = \i] 2 (X),0) and (0 f \S 2 (X) <g> (l f \S 2 (X)w = 
|0, r] 2 (X)). Suppose that the detector outputs are connected so that a single photon produced by either scheme is 
counted. This measurement is described by the projector 1 1/0/) (1/0/ 1 + |0/l/)(0/l/|. The resulting mixed state is 
represented by a density matrix of the form 

P (x) = (1/2) [\ m (x), 0) ( m (x),o\ + |o, m {x)) (o, m (x)\] . (36) 

Its separability is an open question, because a necessary and sufficient condition is known only for mixed systems 
of dimension 2 x 2,2 x 3. However, if we assume that N = 1, i.e., consider a combination of four atoms instead of 
ensembles, then r]2(X) = ^+ = (l/v / 2)(|01) + 1 10) ) and density matrix l|3*H|l describes a four-particle state: 

p(4) = (l/2)(|*+00>(* + 00| + 100*+) (00*+ 1). (37) 

Taking the state of the pair of atoms in the first scheme defined by the two-particle reduced density matrix p(2) = 
(1/2) |\E' + )(\E' + | + 1 00) (00 1 , we can apply the separability criterion. The density matrix partially transposed over the 

variables of the one atom has four eigenvalues one of which is negative 1/4, 1/4, (1 ± v / 2)/4. Therefore, the density 
matrix p(4) is inseparable. 



IX. HIERARCHIC STRUCTURE OF STATES 



Note that expression l|34(l is hierarchically structured. To illustrate this property, we consider a combination of 
schemes generating states of this type. As distinct to schemes using correlation of photocounts, we consider optically 
connected schemes. If an elementary scheme that performs the transformation S n (X) = U n fZ n (W) with X = W (see 
Fig. la) records single-photon output, then the resulting state has the form of (|34H : 

(1/|S„(A> = Vn (X) = t[\XQ . . .0) + • • • + r>0 ... X), (38) 

where w — 1/ y/Prob(l). We define the optical output port of the scheme S n (X) as the one optically coupled to every 
input port. In Fig. la, it is output port 1. The input port of the scheme S n (X) is defined as the optical input port 
of the system of beamsplitters. Then, we can take, for example, p independent schemes represented as (S n (X)) p and 
use their optical outputs as the input of the scheme S p . As a result, we have a new scheme S p ((S n (X)) p ). If it records 
single-photon output, we have an entangled state that consists of lower level entangled states: 

(l f \S p ((S n (X))P)w = VpMX)) = tiMA), . . . 0) + • • • + t p \0, 0, . . . rjn(X)). (39) 

By virtue l|38(l it takes the forme: 

V P (Vn(X)) = Vpn (X). (40) 

Thus, we can formulate the following property. The state r] n (X) defined by l|35(l with n — n\n2-..n p can be 
represented as 

7] n (X) = r) ni (r] n2 {. . . (r] np )). (41) 



This implies that the vector rj n {X) has the structure of an entangled state with respect to any group of s particles, 
where s is such that n/s is a natural number greater than unity. 

When the wavefunction r] n (X) is symmetric, a hierarchically structured representation can be obtained by usin g th e 
permanent expansion defined as a determinant with a summation rule for permutations depending on symmetry 30f . 
In particular, successive decomposition of a determinant with respect to rows or columns and subsequent association of 
summands can be used to represent a permanent in terms of pcrmanents of lower dimension, which reflects hierarchical 
structure. 

For example, when n = 6, it holds that 

V6 (x) = m { m (x)) = m{m(x))). (42) 

This state has the structure of an EPR pair or a W state: 

% (772 (X)) = W(EPR) = EPR{W). 

This example demonstrates that the same state exhibits structure characteristic of entangled states of two different 
types. This property can be used in different applications: the EPR pair can serve as a quantum channel for 
tcleportation or dense coding, while the symmetric W state can be used for cloning. 

To choose a particular structure defined by the dimension of the Hilbcrt space of its element, appropriate basis 
vectors and observables should be introduced. In physical terms, this is equivalent to a two-level approximation. 
Indeed, any group of s particles, where s is such that n/s is a natural number greater than unity, is represented in 
r) n (X) by two states, |0) = S n s (X) = l s . The group can be treated as a two-level particle (qubit) with basis vectors 
S and l s . Such qubits and hierarchically structured states rj n (X) can be used in quantum information processing. 
By analogy with the vector rj n (X) represents a Dicke state only if rj n {X). 



X. CONCLUSIONS 



A model describing resonant interaction of identical two-level atoms with a narrow-band radiation field is used to 
analyze multiparticle entanglement. The interaction is described by an effective Hamiltonian that allows for various 
multiphoton processes. The statistics of radiation and atoms are characterized by a density matrix, for which solutions 
are calculated in secondorder perturbation theory in the interaction strength and exact solutions are found in the 
case of negligible decay. It is shown that the state of any pair of atoms interacting with weak single-mode squeezed 
light is inseparable and robust against decay. It is demonstrated that symmetric entangled multiparticle states can 
be generated by using optical schemes based on singlephoton projection. An optical scheme is described that can be 
used to prepare highly states of entangled atomic ensembles from weakly entangled states by projective measurement. 

This work was supported, in part, by Delzell Foundation, Inc. 
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